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ABSTRACT 

We report on four epochs of observations of the quasar PG 1211+143 using NuSTAR. The net 
exposure time is 300 ks. Prior work on this source found suggestive evidence of an “ultra-fast outflow” 

(or, UFO) in the Fe K band, with a velocity of approximately 0.1c. The putative flow would carry 
away a high mass flux and kinetic power, with broad implications for feedback and black hole-galaxy 
co-evolution. NuSTAR detects PG 1211+143 out to 30 keV, meaning that the continuum is well- 
defined both through and above the Fe K band. A characteristic relativistic disk reflection spectrum 
is clearly revealed, via a broad Fe K emission line and Compton back-scattering curvature. The data 
offer only weak constraints on the spin of the black hole. A careful search for UFO’s show no significant 
absorption feature above 90% confidence. The limits are particularly tight when relativistic reflection 
is included. We discuss the statistics and the implications of these results in terms of connections 
between accretion onto quasars, Seyferts, and stellar-mass black holes, and feedback into their host 
environments. 

Subject headings: galaxies: active - accretion disks - black hole physics - X-rays: binaries - galaxies: 
individual (PG 1211+143) 


1. INTRODUCTION 

The observation of a relation between the masses of 
supermassive black holes at the centers of gal axies and 
the stellar velocity dispersion (M — a re lation; 


Ferrarese 


and Merritt 1 2000 Giiltekin et al.|[2 009) suggests a direct 


link between black holes and their host galaxies. Energy 
and momentum driven out from the central regions push 
gas and dust away, halting star formation and stopping 
AGN fueling ([Silk an d R ees|1998[ |Di M atteo et al.|2005; 
Churazov et al.|2005| see Fabian 2012|±or a review). 'The 
action of AGN feedback could be achieved through the 
powerful radio je ts in the kinetic mode (e.g. McNamara 
and Nulsen|2007|). In the radiative mode, accretion disks 
drive powerful winds that could contribute significantly 
to the energy budget of the BH-galaxy system. 

Observing the properties of such a wind is of great im¬ 
portance, particularly in X-rays where most of the radi¬ 
ation from the expelled material is produced. Although 
warm ab sorber winds are common in the X-ray spectra 
of AGN (|Reynolds||1997| iCrenshaw et al. 1120031 iBlustin 
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et al. 2005), with outflow velocities of ~ 1000 km s 1 
and column densities of log(TVh) ^ 20 — 23 cm -2 , they 
are weak, providing only_^ 0.01% of the AGN bolomet- 


ric luminosity (Blustin et al. 2005| ). The more powerful 


winds seen in several objects with outflow velocities of 
~ 0.1c and column densities of log(TVh) = 24 cm -2 could 


minosity ([Pounds et al. 2003 

|Pounds and Reeves||2009 

Reeves et al 2003; Gappi et a 

. 2009; Ghartas et al. 2009 

Tombesi et al. 

to 

o 

i— 

Ol 

, 2012 Goffordetal. 2013). 


at lea st 35% of observed AGN in X-rays (Tombesi et al. 
2010). However, this number could be an overestimate 
when alternative modeling and more conservative sta¬ 


tistical analyses are c onsidered (Kaspi and Behar 2006 
Vaughan and Uttley| 2008). Establishing how common 
these outflows are, their physical and geometrical prop¬ 
erties is therefore crucial to understanding their contri¬ 
bution to the energy and momentum budget of black 
holes and their hosts. In this le tter, we present analysis 
of the NuSTAR (|Harrison et al.p 

iTAR 1 


2013 


ban* 


observation of the 
(3 — 79 keV) with 


quasar PG 1211+143. NuSTA 
the unprecedented sensitivity at hard (> 10 keV) X-rays 
fixes the continuum and thus allows a meaningful search 
for blue-shifted absorption below 10 keV. 

PG 1211+143 (z = 0.0809) is the archetypical case 
for the ultra-fast outflows in active galaxies. The first 
observation with XMM-Newton in 2001 showed evidence 
for highly blue-shifted absorption lines that are reminis 
nildly relativistic disk winds (^ 0.1c 


cent 


Pounds 


et al. 

2003) 

. The same dataset was analyzed by Kaspi 

and Behar 

(2006) who find a best tit outflow velocity of 

3000 km s _ 

1 instead of the high 24,000 km s -1 . A Chan- 


dra LETG observation showed two redshifted (instead of 
blueshifted!) a bsorption lines at 4.56 and 5.33 keV in the 
source frame ( Reeves et al.|2 005), which, when identified 
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Figure 1 . Left: 3-30 keV light curves from all the four NuSTAR observations of PG 1211+143. The net exposure is shown in each panel. 
Right: The spectrum from the combined observations (combining FPMA and FPMB from all observations) unfolded to a constant to 
remove the effective area of the detector, and plotted as EF(E) using plot eeuf spec in XSPEC. 


as the H-like Ka lines corresponds to inflowing velocities 
of (0.2 — 0.4)c. Later XMM-Newton observations in 2004 
and 2007 showed weaker lines bu t seem to be consist ent 
with the original observations (Pounds and Re eves|2 009), 
or possibly with no absorption lines at all (Tombesi et al 


2010 ) 


2. OBSERVATIONS & DATA REDUCTION 

NuSTAR observed PG 1211+143 in four exposures be¬ 
tween February and July 2014 (The exact dates are: 18 
February 2014, 08 and 09 April 2014 and 07 July 2014). 
The four observations had net exposures 111, 48, 64 and 
74 ks, totaling to nearly 300 ks. The data were reduced 
using HEASOFT V6.16 with the latest calibration (ver¬ 
sion 20141020). We used the scripts nupipeline and 
nuproducts to extract the spectral products. Source and 
background spectra were extracted from regions on (with 
a radius of 2 arcmin) and off source respectively and 
grouped so that there are least 100 source counts per 
bin. The spectra were analyzed using XSPEC V12.8.2. 
Spectral analysis was performed on individual and com¬ 
bined spectra as discussed in section [3j Spectra from the 
two focal point modules A and B (FPMA and FPMB) 
and from different epochs were combined using addspec 
tool in HEASOFT. The response files were combined using 
addrmf with the proper weighting. 

The resulting 3-30 keV light curves from the four ex¬ 
posures are shown in Fig. [l]4eft. The average 3-10 flux 
(from a power-law fit) is 2.7 x 10 -12 ergs cm -2 s -1 which 
is about the same as the first XMM-Newton observation 
of 2001, which had a 3-10 keV flux of 2.5 x 10 -12 ergs 
cm -2 s -1 . The source showed some flux variability be¬ 
tween observations. No strong spectral changes are seen 
apart from a normalization change in the main power-law 
continuum (see section [3|. 

3. SPECTRAL ANALYSIS 

One of the goals of the NuSTAR observation was to 
search for absorption lines from high velocity outflows. 
The spectrum from the new datasets is shown in Fig. [l]- 
right. It has a clear iron K emission line and an excess 
above 10 keV that is most likely due to the Compton 
reflection hump. To be systematic in the search, we con¬ 
sider several baseline models, including simple fits to the 


3-10 keV band so we can directl y compare with the base¬ 
line model used in Tombesi et al. (2010). 


3.1. Searching for Emission and Absorption Features 

In the following discussions, we search for absorption 
(and emission) features by adding a narrovpj Gaussian 
line and doing a systematic scan of the residuals for each 
baseline model. We use a grid of energies between 3.5 
and 9.5 keV in 100 eV steps, and normalization values 
between —lx 10 -5 and 1 x 10 5 (to account for both emis¬ 
sion and absorption). We use a Monte Carlo method to 
obtain signif i cance estimates. As pointed out in |Pro- 


tassov et al.| ( 2002| ), a simple F-test done by comparing 


the improvement in y 2 after adding a Gaussian line is 
not appropriate for obtaining meaningful statistical con¬ 
fidence intervals. The baseline model is not known a 
priori, and neither is the reference (null) distribution of 
the F-statistic. 

This inapplicability of the F-test is not related to the 
quality of data, but rather it is fundamental to any case of 
an added spectral component. Using Monte Carlo meth¬ 
ods provides a way of estimating the unknown reference 
distribution and then using it to assess the significance 
of deviations from the null model. The following steps 
are used to obtain a reference distribution that takes into 
account both the uncertainty in the base line mo del pa¬ 
ramet ers. These are similar to t hose in |Porquet et al. 
(|2004|) and |Tombesi et al.| (|2010|) but additionally take 
into account the baseline model uncertainties: 

A) Start with a baseline model with its best fit Xo an d 
covariance matrix that encodes the uncertainty of best 
value parameters. 

B) Draw N random sets of parameters using the co- 
variance matrix ( N = 1000 in our case). This can be 
achieved for instance by running a Markov Chain Monte 
Carlo (MCMC) parameter search and taking parameters 
directly from the resulting chains (after they have con¬ 
verged). 

C) For each parameter set, fake spectra using the ob¬ 
served response and background files. This produces N 
simulated spectra that are drawn from the baseline model 


10 If a broad component is present in the residuals, several neigh¬ 
boring narrow lines will provide a % 2 improvements. 
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Figure 2. A line search significance plot for five baseline models. The regions with darkest colors indicate that the null model (baseline 
without an additional emission or absorption feature) can be rejected with that significance. The color scale corresponds to log{ 1 — 
significance) and the corresponding significance values are shown in the color bar. The energy axis corresponds to the source reset frame 
energy (z = 0.0809). 


taking into account its uncertainty. 

D) Add a narrow Gaussian line and scan in energy and 
normalization and record the maximum improvement in 
Ay 1 2 . This gives N values of Ay 2 that are used to con¬ 
struct a reference distribution to which observed Ay 2 
values are compared. The significance of an observed 
A xlbs is l-N {Ax 2 >Ax 2 bs) /N, where JV (Ax!>Ax 2 j is the 

number of simulated Ay 2 values greater than A Xobs- 
Several baseline models are used to assess their effect 
on estimating absorption features: 


1. power law + zgauss. This consists of a power-law 
plus a broad Gaussian fitted to the 3-10 keV spec¬ 
tru m similar to those use d in Tombesi et al. (2010) 
and Gofford et al. (2013). 


2. powerlaw+reflionx (Ross and Fabian 2005). 

3. powerlaw+xillver ( Garcia et al.||2014 ). 

4. relxill ( Danser et al.||2013 ). 


For models in 2,3 and 4, the whole NuSTAR band where 
the source is significantly detected above background 
is used. In all models, Galactic absorption is fixed at 
2.7 x 10 20 cm -2 . For the relativistic reflection models, 
we assume a power-law emissivity and we fit for its in¬ 
dex ( q ) along with the black hole spin, photon index, 
ionization parameter, iron abundance, inclination, high 
energy cutoff and reflection fraction. The details of the 
model parameters are shown in Table [l] For each base¬ 
line model, we repeat steps A-D. The results are shown 
in Fig. [2] 


3.2. Line Search Results 

It is clear from Fig. [2] that there is no significant ab¬ 
sorption feature that is persistent with all baseline mod¬ 
els. The phenomenological description of the first model 
gives an excellent description to the data in the 3-10 
keV band with no significant residuals. In the ref lionx 
fit, an absorption feature at 7.1 keV is apparent along 
with the emission feature at 5.8 keV. The same feature at 
5.8 keV persists for the xillver fit. Both the ref lionx 
and xillver models account for the Fe Ka emission 



powerlaw+zgauss (3-10 keV) 

T'zgauss 

6.28 ± 0.07 keV 

^zgauss 

0.33 ± 0.07 keV 


powerlaw+reflionx (3-50 keV) 

Ap e (solar) 

0.7 ±0.1 

r 

2.45 ± 0.06 

logit) 

1.9 ±0.2 


powerlaw+xillver (3-50 keV) 

Ap e (solar) 

0.5 ±0.1 

r 

2.50 ±0.06 

i°g(£) 

1.3 ±0.4 


relxill (3-50 keV) 

q 

2.2 ±.5 

a 

-0.13 ±0.8 

incl. 

28 ±7 

r 

2.51 ±0.2 

log(0 

1.3 ±0.9 

Ap e (solar) 

0.7 ±0.1 

£cut 

124 (keV) (lower limit) 

Ref. fract. 

9 Q+ 1 - 6 


Table 1 

Fit parameters for the baseline models discussed in section |±T] 


at 6.4 keV. The residuals at ^ 5.8 indicate that addi¬ 
tional emission is required red-ward of the line, suggest¬ 
ing line broadening. If broadening is included in the fit 
(relxill, relxill+xillver), the feature disappears. 

The absorption in the ref lionx fit at 7.1 keV disap¬ 
pears too when relativistic blurring is included. The last 
model (relxill+xillver) is included for completeness. 
It models reflection from both the inner (with relativis¬ 
tic broadening) and outer regions. Although, the addi¬ 
tional component improves the y 2 slightly compared to 
relxill alone, its requirement is not very significant. 
The conclusion from Fig. [2] is that absorption features 
depend on the chosen continuum, and when the contin¬ 
uum is pinned down with NuSTAR coverage at high en¬ 
ergies and the inclusion of relativistic reflection, no sig¬ 
nificant absorption is seen. 

The analysis summarized in Fig. [2] was done using a 
combination of all the spectra (FPMA and FPMB from 
four observations) to obtain the best signal. We tested 
for any possible flux dependency by grouping observa- 
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tions 1 and 4 together and 2 and 3 together (see light 
curve in Fig. [I]4eft), and the results do not change; no 
additional absorption lines are required by the data. 


3.3. Comparison with XMM-Newton 

The first XMM-Newton observation in 2001 showed an 
absorption line at ~ 7.1 keV in the observe r’s frame p lus 
possi ble additional lines at lower energies (Pounds et al. 
2003). We reanalyzed this XMM-New ton d ataset using 
the same method discussed in section 13.11 and find the 
line at ^ 7.5 keV to be significant at the ~ 99% level 
with an equivalent width of 104 =b 52 eV. The flux levels 
of the first XMM-Newton observation and the NuSTAR 
observation discussed here are comparable; therefore, we 
also assess whether an absorption line similar to that 
seen in XMM-Newton would have been detectable in the 
current NuSTAR data. We fitted the XMM-Newton PN 
spectrum with a model consisting of a power-law, rela¬ 
tivistic reflection and an absorption line. We then used 
this model to f ake N uSTAR data, and followed the proce¬ 
dure of section [3T] to search for absorption line. The line 
is detected significantly in this simulated data showing 
that an absorption line would have been easily (> 8a sig¬ 
nificance) detected if it was present in the data. Adding 
a Gaussian line to our best fit model (relxill) with the 
energy and the width of absorption line seen in the XMM- 
Newton spectrum provides a fit improvement of Ay 2 = 3 
for one degree of freedom. The equivalent width of this 
line is 22 d= 22 eV. 

Subsequent to the first observation, XMM-Newton ob¬ 
served PG 1211+143 several times, and when we analyze 
these observations self-consistently, none of them show 
any significant absorption features when using a baseline 
model consisting of a power-law and a relativistic line. 


3.4. Relativistic Reflection Parameters 

The best fit model to the NuSTAR data, shown in 
Fig. [3J consists of an absorbed power-law and relativis¬ 
tic renection, both modeled with relxill. Only data be¬ 
low 10 keV are shown to highlight the differences among 
models. As was clear from Fig. [2j the data show an ex¬ 
cess at ^ 5 keV when fitted with simple reflection models. 
The addition of relativistic blurring (relxill) accounts 
for the excess. The data, however, do not allow very 
strong constraints on the spin of the black hole. The 
parameters of the best fit model are shown in Table [l] 

If we fit for the inner radius of emission instead of the 
black hole spin, we obtain the marginal probability con¬ 
tours shown in the bottom panel of Fig. [3j This figure is 
obtained by running a Monte Carlo Markov Chain on 
the best fit relxill to obtain posterior probability dis¬ 
tributions of the fit parameters and then marginalize over 
the parameters that are not of interest. The inner radius 
of emission and emissivity index are degenerate in the fit, 
but it is clear that although the extra broadening of the 
Fe Ka line is required, it is not very large. The inner ra¬ 
dius of emission ranges from ^ 4 to 40 gravitational radii 
(r g = GM/c 2 ) depending on the disk emissivity. The so¬ 
lution with a steep emissivity and large radii might be 
unphysical, and although this data do not allow us to 

11 We use xspec_emcee available from: 
https: / / github .com / j eremysanders / xspec.emcee 
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Figure 3. Top: The spectrum of PG 1211+143 along with the 
best fit baseline models discussed in section |TT] plotted in the 3-10 
keV only where the models differ. The x-axis represents the rest- 
frame energy of the source. Bottom: Marginal probability density 
contours in the relxill model for the inner radius of emission and 
the emissivity index. The contours are for the 68, 90 and 99% 
confidence levels. 


firmly locate the X-ray emission, the detection of rever¬ 
beration time delays in this objects in the XMM-Newton 
data (de Marco et al. 2011) would favor the solution with 
a small radius for the emission. 


4. DISCUSSION AND CONCLUSIONS 

We have analyzed a deep 300 ks observation of the 
quasar PG 1211+143, obtained with NuSTAR in four 
time intervals spanning six months. When fit with a 
simple power-law model, the total spectrum shows clear 
signatures of disk reflection (Fig. fright). Fits with the 
newest, most sophisticated reflection models currently 
available provide significantly improved fits. The resid¬ 
uals to fits with simple reflection show an excess around 
~ 5 keV indicating that relativistic blurring is required, 
clearly signaling that X-ray emission probes the very in¬ 
nermost portion of the accretion flow in PG 1211+143. 
The data do not permit a strong constraint on the spin 
of the black hole, but when reverberation delays seen 
previously in XMM-Newton data are considered, a solu¬ 
tion with a small inner radius (~ 4 r g , corresponding to 
a black hole spin of ~ 0.5) is supported by the data. We 
thoroughly examined the NuSTAR data both in segments 
and in total, and the spectra show no compelling evi- 
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dence of a steady or variable ultra-fast outflow seen pre¬ 
viously in this source. The limits are particularly tight 
once disk reflection has been modeled. 

NuSTAR clearly has the sensitivity and spectral res¬ 
olution required to detect extraordinary outflows, when 
they are present. Observations of PDS 456, for instance, 
reveal a likely P-Cygni profile in that quasar, signaling a 
massive and powerful outflow from that black hole into 
its host galaxy (Nardini et al. 2014, submitted). NuS¬ 
TAR spectra from observations of NGC 1365 were able 
to disentangle variable absorption, outflows, and disk re¬ 
flection in that source, leading to a robust measure of the 


black hole spin (Risaliti et al. 2013 Walton et al. 2014). 
At the opposite end of the mass scale, NuSTAR observa- 
tions of the black hole candidate 4U 1630—472 revealed 
strong absorption that ma y be a UFO -like outflow, in 
addition to disk reflection (King et al. 2014). The true 
extent of these ultra-fast outflows and their persistence 
is unclear. 

The non-detection of ultra-fast outflows in the proto¬ 
typical object PG 1211+143 in the current observation 
(and also previous observations subsequent to the first 
XMM-Newton detection) raises several questions on how 
common these might be in this and other sources. Sys¬ 
tematic searches suggest that ^ 40% of objects o bser ved 


with XMM - Newton and Suzaku show them (Tombesi 


et al. 2010 Gofford et al. 2013). This number could 


be much lower if a careful statistical search is performed 
that takes into account physical baseline models along 
with their observational uncertainties. This comes in the 
form of including the best fit uncertainties in simulating 
spectra during the Monte Carlo tests and also allowing 
for that uncertainty when fitting for narrow features at 
arbitrary energies in the spectrum. 

The true extent of UFO’s could be affected by the 
significance of any single detection, which we have dis¬ 
cussed, and also by their variability. A non-detection 
could simply be due to the fact they are tra nsient in 
nature as sugg ested by the the simulations ( Proga and] 


Kallman 2004). The case of PG 1211+143 shows that 
the UFO is seen in one out of eleven observations ( XMM- 
Newton , Suzaku , Chandra and NuSTAR ) between 2001 
and 2014. With the caveat that we still have small num¬ 
bers, if one assumes the observations are done randomly 
as far as the AGN is concerned, it appears that a variabil¬ 
ity argument would also suggest that the contribution of 
UFO’s is an order of magnitude lower than previously 
inferred. 
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